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1.  Introduction 


Over  the  past  decade,  there  has  been  considerable  scientific  interest  focused  around  the  unique  properties 
and  device  applications  of  single  sheet  and  multilayered  van  der  Waals  materials  [1-6].  When  exfoliated  or 
grown  as  single  atomic  layers,  many  of  these  two-dimensional  (2D)  materials  exhibit  electronic  and 
optoelectronic  properties  which  differ  fundamentally  from  those  of  the  parent  bulk  materials  [7,  8].  These 
changes  are  indicative  of  the  small  but  appreciable  electronic  coupling  that  arises  between  the  neighboring 
layers  of  these  bulk  crystals.  For  instance,  single -layered  graphene  exhibits  the  quantum  Hall  effect  at  room 
temperature.  This  is  not  seen  in  its  bulk  counterpart,  graphite,  and  demonstrates  graphene’s  high-mobility 
two-dimensional  electronic  structure  [9].  Additionally,  the  normally  indirect  band  gaps  of  the  group  V1B 
layered  transition  metal  dichalcogenides,  such  as  M0S2,  transform  into  direct  band  gaps  only  when 
exfoliated  or  grown  as  single  layers  [10-13].  This  dramatic  contrast  between  single-  and  multi-layered  2D 
crystals  implies  that  there  exists  significant  coupling  between  the  different  layers.  In  contrast,  a  material 
where  this  interlayer  coupling  is  strongly  suppressed  or  can  be  manipulated  would  enable  observation  of 
these  exotic  2D  phenomena  even  in  relatively  thick  multilayered  crystals. 

The  Si/Ge/Sn  graphane  analogues  are  a  unique  class  of  van  der  Waals  materials  that  feature  a  2D 
electronically  active  main  group  lattice  in  which  every  atom  is  terminated  with  a  covalent  ligand  [14-17]. 
This  covalent  termination  serves  to  dramatically  suppress  band  dispersion  in  the  out-of-plane  direction  as 
compared  to  previously  studied  van  der  Waals  materials.  Furthermore,  varying  the  identity  of  the  covalently 
terminated  ligand  enables  tuning  of  the  in-plane  and  out-of -plane  electronic  and  thermal  transport  behavior, 
optical  properties,  thermal  stability,  and  topological  phases  in  these  materials  [16,  18-20].  The  recent 
development  of  a  facile  approach  to  the  synthesis  of  Ge  graphane  analogues  has  resulted  in  mm-scale 
crystals  of  both  hydrogen-terminated  germanane  (GeH)  and  methyl -terminated  germanane  (GeCH3)  [14, 
16,  21],  Hydrogen-terminated  germanane  is  of  particular  interest  due  to  the  predictions  of  a  room 
temperature  phonon -limited  electron  mobility  of  18,000  cm2/V-s — a  value  approximately  five  times  larger 
than  the  mobility  in  bulk  Ge — and  a  direct  bandgap  of  1.59  eV.  Density  functional  theory  predicts  little 
change  in  the  electronic  structure  between  an  isolated  sheet  of  GeH,  and  the  two-layer  unit  cell  [16]. 
Consequently,  there  is  only  a  weak  electronic  coupling  between  crystal  layers.  Therefore,  the  high  mobility 
and  direct  band  gap  are  predicted  to  exist  both  in  exfoliated  flakes  and  in  the  bulk  crystalline  material. 
While  studies  on  electronic  transport  in  silicene  and  related  silicon  based  graphane  analogues  exist,  to  date 
there  have  been  none  on  the  germanium  based  graphane  analogues  [22,  23]. 

Herein,  we  measure  for  the  first  time  the  electrical  transport  properties  of  GeH  through  measurements 
on  both  undoped  and  phosphorous  doped  GeH  multilayer  crystals.  Dopants  are  incorporated  into  the  GeH 
lattice  by  first  integrating  stoichiometric  equivalents  of  the  elemental  dopants  into  the  precursor  CaGe2  Zintl 
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phase.  The  effectiveness  of  this  phosphorous  doping  is  verified  by  observation  of  a  reduction  in  sheet 
resistance  by  up  to  three  orders  of  magnitude.  There  is  also  a  significant  anisotropy  between  electronic 
transport  parallel  and  perpendicular  to  the  layered  structure.  Additionally,  there  is  a  strong  atmospheric 
dependence  in  the  resistance,  with  the  sample  resistance  decreasing  when  water  vapor  is  present  in  the 
environment.  These  results,  combined  with  a  finite  element  analysis  of  the  transport  in  GeH,  are  used  to 
extract  the  in-plane  and  out-of-plane  resistivities  and  to  model  the  current  distribution  in  GeH  crystals.  This 
suggest  that  the  sample  current  flows  primarily  through  the  outermost  layers  due  to  dopant  activation  by 
the  water  vapor  present  in  the  environment.  This  dopant  activation  results  in  a  modeled  resistivity  four 
orders  of  magnitude  smaller  than  those  calculated  in  vacuum. 

2.  Experimental  Methods 

Single-crystalline  flakes  of  GeH  were  synthesized  using  methods  previously  reported  [8].  In  brief,  the 
layered  CaGea  Zintl  phase  precursor  was  produced  by  annealing  a  stoichiometric  ratio  of  calcium  and 
germanium  in  a  sealed  quartz  tube  at  950-1050  °C  under  vacuum.  The  samples  were  subsequently  cooled 
over  the  course  of  2-10  days.  The  CaGe2  crystals  were  then  converted  to  GeH  by  treatment  in  concentrated 
HC1  for  5-14  days  at  -40  °C.  This  process  yielded  GeH  crystals  which  measure  up  to  4  mm  laterally  and 
have  a  thickness  of  5-20  pm.  To  tune  the  electronic  properties  of  the  resulting  GeH,  phosphorous  was  used 
as  an  extrinsic  n-type  dopant.  Phosphorous  doped  GeH  crystals  (P:GeH)  were  obtained  by  incorporating 
stoichiometric  amounts  of  phosphorous — 0.5  and  1.0%  atomic  percentage — to  the  initial  Ca/Ge  mixture, 
followed  by  the  same  annealing  and  HC1  treatment  processes  used  for  the  undoped  GeH  synthesis. 

The  crystallinity  of  the  resulting  GeH  and  P:GeH  crystals  was  confirmed  via  powdered  x-ray  diffraction 
(XRD),  Raman  spectroscopy,  and  scanning  electron  microscopy  (SEM).  XRD  was  performed  on  a  Bruker 
D8  powder  X-ray  diffractometer  with  Cu  Ka  radiation.  The  Raman  scattering  spectra  were  collected  on 
Renishaw  InVia  Raman  equipped  with  a  CCD  detector  using  a  633  nm  He-Ne  laser  for  excitation.  The 
SEM  micrographs  were  taken  using  a  FEI  Helios  Nanolab  600  dual  beam  focused  ion  beam  and  scanning 
election  microscope  at  current  of  43  pA  and  voltage  of  10  kV. 

Top-contacted  GeH  samples  were  fabricated  by  electron-beam  deposition  of  Ag/Au  (100/20  nm) 
through  a  shadow  mask,  creating  a  25  pm  channel  between  the  two  contacts.  Depositions  were  conducted 
with  the  sample  at  room  temperature  to  avoid  amorphization  which  occurs  above  75  °C  [16].  Silver  (with 
a  Au  capping  layer  to  prevent  oxidation)  was  selected  as  a  contact  metal  after  showing  that  the  silver 
contacts  on  P:GeH  had  smaller  resistances  and  linear  I-V  curves  indicating  a  low  contact  resistance  and 
minimal  Schottky  barrier.  Other  metals  tested  were  Ti/Au,  Au,  Pd,  Cr,  and  Pt.  Connections  to  the  electron 
beam  deposited  sample  contacts  were  made  by  connecting  thin  Au  wires  to  the  deposited  contacts  using  Ag 
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epoxy.  The  out-of-plane  samples  were  fabricated 
by  placing  a  drop  of  Ag  epoxy  on  a  Si  02/Si 
substrate  (300nm  SiOi),  setting  the  GeH  flake  on 
top  of  this  epoxy,  and  then  placing  a  second  epoxy 
spot  on  the  top  of  the  flake.  Here  again,  Au  wire 
was  used  to  make  contact  to  the  epoxy  contacts. 

Transport  current-voltage  scans  were 
measured  using  a  Keithley  4200-SCS,  using  a 
probe  station  manufactured  by  Lake  Shore 
Cryotronics  Inc.  (Westerville,  OH).  Unless 
otherwise  noted,  all  I-V  curves  are  measured  in 
ambient  atmosphere,  at  room  temperature,  and  in 
a  dark  environment  (to  prevent  photo-excited 
carrier  contributions  to  the  conductivity). 

3.  Results  and  Discussion 

3.1  Sample  Characterization 

The  purity  of  the  resulting  P:GeH  crystals  is  confirmed  through  X-ray  crystal  diffraction  (XRD)  and  Raman 
spectroscopy  (figure  1(a)  and  (b)).  The  samples  can  be  indexed  to  a  hexagonal  unit  cell  with  a  =  3.880  A 
and  c  -  1 1.04  A  (two  GeH  layers  per  c-unit  cell  spacing),  consistent  with  previous  measurements  [16,  21]. 
The  XRD  spectra  indicate  that  there  are  no  additional  impurity  phases  or  significant  changes  in  the  lattice 
parameters  for  the  doped  samples.  The  Raman  spectra  of  the  GeH  and  P:GeH  samples  show  the  expected 
E2  mode  at  302  cm  1  as  well  as  the  Ai  mode  at  228  cm'1  with  no  discernable  changes  in  peak-shape  or 
position  [16,  21].  Scanning  electron  micrographs  (figure  1(c),  (d),  and  (e))  of  the  GeH  crystals  confirm  that 
they  are  smoothly  faceted  with  no  visible  gaps  between  the  GeH  sheets. 

3.2  Atmospheric  Dependence  and  Anisotropy  of  Resistance  Measurements 

In  bulk  crystalline  samples  of  2D  materials,  there  is  typically  a  strong  anisotropy  (ranging  from  2-6  orders 
of  magnitude)  between  resistivity  in  the  plane  of  the  2D  sheets  and  between  the  sheets  [23-26].  In  order  to 
probe  this  anisotropy,  as  well  as  the  effectiveness  of  the  P  dopants,  two  different  contact  geometries  were 
employed  to  measure  current-voltage  (I-V)  scans.  The  in-plane  (IP)  conductivity  was  probed  through  the 
use  of  a  top  contact  (TC)  geometry  (figure  2(a),  2(b)  inset),  while  an  out-of-plane  (OOP)  geometry  (figure 
2(c),  2(d)  inset)  was  used  to  measure  the  interplane  resistivity. 


Figure  1.  Raman  (a)  and  XRD  (b)  spectra  of  HC1 
deintercalated  GeH  with  0.0%  (undoped,  black),  0.5% 
(red),  and  1.0%  (blue)  phosphorous  incorporation  in  the 
C’aGe,  precursor,  (c-e)  SEM  images  of  the  bulk  GeH 
crystals. 
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Representative  I-V  curves  for  TC  and  OOP 
geometries  are  shown  in  figures  2(b)  and  (d), 
respectively  (data  shown  here  for  0.5%  P:GeH). 
In  addition  to  the  large  difference  in  resistances 
between  the  two  contact  geometries  (extracted 
from  the  slopes  of  the  IV  curves),  there  is  also  a 
variation  in  the  sensitivity  of  the  two 
measurements  to  the  environmental  conditions.  In 
particular,  measurements  conducted  on  TC 
devices  show  an  increase  in  current  by  three 
orders  of  magnitude  when  exposed  to  the  ambient 
atmosphere  compared  to  measurements  under  a 
vacuum  of  5  x  10'4  mbar.  In  contrast,  the  OOP 
devices  experience  almost  no  change  in  resistivity 
as  both  vacuum  and  ambient  atmosphere 


Figure  2.  Schematics  of  the  top-contact  (a)  and  out-of- 
plane  (c)  geometries  used  to  study  the  transport  properties 
in  GeH  crystals,  (b)  and  (d)  show  representative  IV  curves 
both  in  ambient  atmospheric  conditions  (blue)  and  in  a 
vacuum  of  5 x  10 4  mbar  (black).  Out-of-plane 
measurements  are  at  the  lower  limit  of  the  hardware  setup. 
Insets  show  characteristic  micrographs  of  samples  (e) 
Photocurrent  as  a  function  of  light  energy  which  matches 
the  absorbance  curve  from  Bianco  et  al.  showing  an 
increase  in  conductivity  with  exposure  to  above  bandgap 
light  [16], 


measurements  exhibit  currents  at  the  lower  limit 
of  the  hardware  setup.  This  lack  of  variation  could 
be  due  to  an  absence  of  atmospheric  sensitivity  in 
the  OOP  conduction,  or  the  magnitude  of  this 
change  is  beyond  our  current  sensitivity.  To 
confirm  that  the  observed  TC  geometry’s  low 
current  signals  are  generated  by  transport  through 
the  P:GeH  crystal  rather  than  through  parallel 


conduction  channels  such  as  the  substrate  or  a  surface  water  layer,  variable  wavelength  photoconductivity 


measurements  were  performed  at  light  energies  of  1.2-2. 2  eV  under  vacuum.  Figure  2(e)  shows  the 
photocurrent  signal  (black  data)  on  top  of  the  diffuse  reflectance  absorption  measured  previously  by  Bianco 


et  al.  in  GeH  (red  data)  [16].  The  increase  in  photocurrent  at  wavelengths  above  the  1.59  eV  bandgap 
correlates  well  with  the  absorption  spectra  and  confirms  that  the  current  measured  is  due  to  conduction 
inside  the  P:GeH  samples. 


Since  the  photoconductivity  measurement  confirms  that  the  conduction  is  occurring  in  the  GeH  crystal, 
the  atmospheric  dependence  observed  in  the  TC  geometry  suggests  that  a  component  of  the  ambient 
atmosphere  may  be  activating  dopants  on  the  surface  of  the  P:GeH  crystals.  Figure  3  shows  the  positive 
voltage  I-V  scans  for  a  prototypical  0.5%  P:GeH  sample  in  the  TC  geometry  upon  exposure  to  a  variety  of 


5 


atmospheres.  Negative  scans  show  similar  behavior 
and  can  be  found  in  the  supplemental  material 
(figure  SI).  The  sample  was  first  measured  in  air 
and  then  the  chamber  was  pumped  down  to  a 
vacuum  of  5  x  10'4  mbar  for  a  second 
measurement.  The  chamber  was  subsequently 
backfilled  to  atmospheric  pressure  with  Ar,  O2,  and 
32  mbar  water  in  argon  mixture  (tbO/Ar).  The 
system  was  pumped  down  to  the  base  pressure  after 
each  gas  exposure.  On  exposure  to  both  air  and 
water  vapor,  the  maximum  current  increased  as 
compared  to  the  vacuum,  Ar,  and  O2  I-V  scans. 
These  measurements  demonstrate  that  the 
component  of  the  ambient  atmosphere  responsible 
for  the  decreased  resistivity  signal  in  the  TC 
measurements  is  the  water  vapor. 

3.3  Doping  Dependence 

Next  the  changes  in  electronic  properties  in  samples  with  different  phosphorous  incorporation  were  probed 
in  both  the  TC  and  OOP  geometries.  We  investigate  three  cases:  undoped  GeH,  0.5%  P  incorporation,  and 
1 .0%  P  incorporation.  Figure  4  shows  the  respective  I-V  curves  at  different  phosphorous  incorporations  for 
the  TC  geometry  (figure  4(a))  and  the  OOP  geometry  (figure  4(b)).  The  TC  devices  show  an  increase  in 
resistance  by  three  orders  of  magnitude  under  increased  phosphorous  doping,  while  the  OOP  geometry 
show  no  apparent  increase.  Here  again,  all  of  the  OOP  measurements  show  currents  that  are  at  the  hardware 
limits.  The  variation  in  sheet  resistance  of  the  different  phosphorous  incorporation  regimes  is  shown 
quantitatively  in  figure  4(c)  for  the  TC  measurement  geometry.  Resistances  were  extracted  from  the  slope 
of  the  best  fit  line  for  each  I-V  with  three  to  five  samples  measured  at  each  P  incorporation  value.  In  order 
to  account  for  sample-to-sample  geometric  effects  in  the  extracted  resistances,  the  resistances  were 
converted  to  sheet  resistances.  Error  bars  represent  the  sample  to  sample  variation  observed  in  the  doping 
series  measurements. 

While  the  large  resistance  anisotropy  of  the  TC  and  OOP  contact  geometries  is  a  useful  metric,  the 
intrinsic  in-plane  and  out-of-plane  resistivities  (px  and  pz,  respectively)  and  the  anisotropy  ( A  =  pz/px) 
between  them  provide  insight  into  the  inter-  and  intra-layer  coupling  of  the  material.  The  OOP  resistivity 


Figure  3.  Representative  positive  voltage  I-V  curves 
showing  the  atmosphericic  dependence  of  P:GeH  using 
top  contact  geometry.  Samples  was  exposed  to  air 
(blue),  water/argon  mixture  (green),  oxygen  (red), 
argon  (violet),  and  vacuum  (black).  The  sample  was 
pumped  down  to  a  base  vacuum  of  5  *  1CT4  mbar 
between  measurements. 
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(pz)  can  be  extracted  directly  from  the  OOP 
measurement  geometry  using  the  standard 
equation  for  uniformly  distributed  current  flow 
inside  a  material: 

(WL\ 

Pz  —  ^oop  J 

Where  Roop  is  the  resistance  measured  in  the 
OOP  geometry,  W  and  L  are  the  two  in-plane 
sample  dimensions,  and  T  is  the  out  of  plane 
thickness  of  the  sample.  Since  the  OOP 
geometry  measurements  exhibit  a  current  signal 
at  the  limit  of  our  hardware,  only  a  lower  bound 
of  pz  >  5xl014  O  •  cm  can  be  calculated  for  all 
three  doping  regimes.  This  value  is  comparable 
to  that  measured  previously  in  half  hydrogen  and 
half  hydroxide  terminated  silicane 
[SiH0.5(OH)o.5]  at  250  °C  [23], 

On  the  other  hand  px  cannot  be  directly 
extracted  from  the  TC  geometry  measurements 
due  to  the  non-uniformity  of  the  current  distribution  which  is  a  combination  of  in-plane  and  out-of-plane 
current  paths.  Further,  the  strong  atmospheric  dependence  observed  in  the  TC  geometry  implies  transport 
in  this  geometry  is  more  sensitive  to  the  surface  of  the  sample  (and  hence  px)  than  it  is  to  the  bulk  of  the 
sample  (and  pz).  This  suggests  an  anisotropy  between  px  and  pz  that  is  partially  due  to  the  atmospheric 
activation  of  dopants  present  in  the  samples.  In  order  to  probe  this  current  distribution  in  the  TC 
measurement  and  to  place  bounds  on  px  (and  hence  the  anisotropy  between  px  and  pz),  we  developed  a 
finite  element  model  of  transport  inside  of  multilayered  2D  materials. 

3.3  Finite  Element  Modelling  of  In-Plane  Resistivity 

The  finite  element  model  used  here  is  illustrated  in  figure  5(a).  The  bulk  crystal  is  modeled  by  a  network 
of  in-plane  and  out-of-plane  resistors  with  resistances  Rx  and  R7.  respectively.  This  simple  model  can  be 
modified  to  incorporate  two  materials  with  different  internal  resistance  values  by  stacking  blocks  of 
materials  on  top  of  each  other,  each  with  their  own  internal  resistances.  We  use  primed  variables  to  refer 
to  the  resistances  ( R'x  and  /?')  in  the  top  material  block,  and  unprimed  variables  (Rx  and  Rz)  to  note  the 


(a)  Top  Contact  (b)  Out-of-Plane  Contact 


(c) 


Figure  4.  Doping  dependence  of  GeH  samples  with  0.0% 
(black  data),  0.5%  (red  data),  and  1.0%  (blue) 
incorporation  of  phosphorous  in  the  CaGe,  precursor. 
Representative  I-V s  in  the  top-contact  (a)  and  out-of-plane 
(b)  geometries.  Out-of-plane  I-Vs  show  a  signal  below  the 
limits  of  the  hardware  setup,  (c)  Sheet  resistances  in  top- 
contact  geometry  samples  for  different  doping  regimes. 
Error  bars  represent  the  sample  to  sample  variation. 
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resistances  in  the  lower  block.  This  network  of  resistors  can  then  be  solved  for  the  equivalent  sheet 
resistance  of  the  complete  network.  The  values  of  network  resistors  are  tuned  in  order  to  match  this 
equivalent  sheet  resistance  to  the  measured  sheet  resistance.  Additionally,  this  model  yields  the  voltage  at 
each  node  and  the  current  across  the  resistors  between  nodes.  The  input  model  resistances  Rx  and  Rz  for 
both  layers  can  be  directly  related  to  the  physical  resistivities  of  the  materials  of  the  different  blocks  ( px 
and  pz)  through  the  relations: 

Rx  =  px  (rd  and  R?  =  pz  (rd 

where  l,  t,  and  w  are  the  physical  cell  dimensions  used  to  describe  the  model.  In  the  case  of  GeH,  we  let 
l  =  t  =  0.6  nm — the  thickness  of  a  physical  GeH  sheet — and  let  w  be  the  width  of  the  samples  themselves. 
This  network  of  resistors  is  solved  analytically  using  Kirchhoff  s  circuit  laws  to  generate  a  set  of  linear 
equations  which  are  solved  for  the  voltage  at  each  node,  the  current  across  each  resistor,  and  the  equivalent 
sheet  resistance  of  the  network  as  a  whole.  The  model  was  implemented  using  National  Instruments  (NI) 
Lab  VIEW  software  along  with  Nl’s  multicore  analysis  and  sparse  matrix  toolkit.  We  apply  this  model  to 
the  0.5%  P:GeH  measurements. 

To  extract  the  in-plane  resistivity  of  the  P:GeH,  the  pz  value  extracted  from  OOP  measurements  is  held 
constant,  while  px  is  iterated  until  the  equivalent  sheet  resistance  of  the  model  matches  the  sheet  resistance 
measured  in  the  TC  geometry  (RT c).  The  vacuum  measurements  (black  curves  in  figures  2(b)  and  2(d))  are 
simply  modeled  by  a  single  block  of  material.  Using  the  sheet  resistance  of  i?TC  =  1.0  x  1 0 14  Q/n  measured 
in  the  TC  geometry  and  pz  =  5  x  1 0 1 4  Q  •  cm,  extracted  from  OOP  geometry  measurement,  the  in-plane 

n 

resistivity  is  calculated  to  be  px  =  1 .4  x  1 0  O  ■  cm.  This  resistivity  is  four  orders  of  magnitude  smaller 
than  the  in-plane  resistivity  measured  in  the  half  hydrogen  and  half  hydroxide  terminated  silicane 
[SiHo.5(OH)o.s]  at  250  °C,  and  twelve  orders  of  magnitude  smaller  than  the  extrapolated  room  temperature 
value  [23]. 

To  model  the  resistivities  of  the  samples  exposed  to  water  vapor  (I-V  measurements  shown  in  blue 
curves  in  figures  2b  and  2d),  the  two  block  model  is  used.  Here,  the  top  block  is  made  up  of  a  number  of 
surface  layers  in  which  the  atmosphere  has  water  activated  carriers  (assigned  values  of  px  and  pz)  while  the 
second  block  is  the  non-activated  bulk  (assigned  values  of  px  and  pz).  The  value  for  px  and  pz  are  the  same 
values  extracted  from  the  vacuum  measurements  (px=  1.4  x  1 07  O  ■  cm,  pz  =  5  x  10 14  Q  •  cm).  The  values 
of  px  and  pz  are  then  varied  so  that  the  equivalent  sheet  resistance  of  the  network  matches  the  observed 
sheet  resistance  of  RTC  =  3 . 8  x  1 0 1 0  Q/n.  The  number  of  water  activated  layers  is  varied  in  the  model  as  a 
means  of  accounting  for  the  surface  roughness  in  the  measured  samples. 


Z  (Layers)  z  (Layers) 


(a) 


Y®— ->x 

j.  K 


Water  Activated  Layers 


■» 

• 

• 

• 

• 

♦ 

A 

Insulating  Bulk  Layers 


X  (pm)  X  (nm)  x  (h™) 

Figure  5.  (a)  Cartoon  showing  the  model  developed  to  investigate  current  flow  in  a  layered,  multi-material  structure. 
Material  layers  are  split  into  cells  with  dimensions  equal  to  that  of  a  single  layer  thickness.  Material  resistances  (in¬ 
plane  and  out-of-plane)  are  input  as  the  pairs  /?',  /?'  and  Rx,  Rz,  and  the  network  of  resistors  is  solved  for  the  voltage 
at  each  voltage  node  (V  )  and  the  equivalent  resistance  of  the  network.  Voltage  and  current  maps — top  and  bottom 
of  (b),  (c),  and  (d),  respectively — are  shown  for  a  smooth  surface  (a),  a  rough  surface  with  anisotropic  conductivity 
in  the  water  activated  0.5%  P:GeH  block  (b),  and  a  rough  surface  with  isotropic  conductivity  in  the  water  activated 
P:GeH  block  (c).  The  smooth  surface  is  simulated  using  one  single  layer  of  water  activated  P:GeH,  while  rough 
surfaces  are  simulated  with  ten  layers  of  water  activated  P:GeH.  Dashed  lines  distinguish  between  the  upper  water 
activated  block,  and  the  bottom  non-water  activated  block  of  the  model. 
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Rougher  samples  are  expected  to  have  a  higher  number  of  layers  exposed  to  the  atmosphere  while  an 
atomically  smooth  sample  would  only  have  the  single  topmost  layer  exposed  to  the  atmosphere.  For  this 
reason,  a  rougher  sample  is  modeled  assuming  the  top  ten  layers  of  GeH  are  water  activated,  while  a  smooth 
sample  assumes  only  the  top  layer  is  water  activated. 

For  the  case  of  an  atomically  smooth  sample  we  assume  that  only  a  single  layer  of  material  experiences 
water  activation  of  dopants,  thus  there  is  no  p'z  value  and  only  p'x  is  varied  until  the  effective  sheet  resistance 
matches  the  observed  sheet  resistance.  This  occurs  at  p'=  2.3  x  1  ()  ’  O  ■  cm.  The  resultant  current  and 
voltage  maps  for  this  case  are  shown  in  figure  5(b).  Modeling  indicates  that  in  this  case,  >99%  of  the  current 
travels  solely  in  this  topmost  layer.  The  modeled  p'x  value  is  nearly  four  orders  of  magnitude  smaller  than 
the  value  derived  in  the  unactivated  case  ( px  =  1.4  x  10  Q  •  cm),  and  eleven  orders  smaller  than  the  out 
of  plane  resistivity  (pz  =  5*1014  G-  cm). 

A  more  realistic  view  of  the  P:GeH  samples  is  that  the  surface  is  not  atomically  flat.  In  this  scenario, 
several  layers  are  exposed  to  the  atmosphere  and  experience  dopant  activation.  To  model  this,  the  top  ten 
layers  of  the  sample  are  assumed  to  be  water  activated.  This  corresponds  to  a  roughness  of  6  nm.  While 
water  is  unlikely  to  diffuse  through  the  bulk  of  the  GeH  samples  and  change  the  OOP  resistance 
measurements,  the  measured  currents  in  the  OOP  geometry  are  below  the  minimum  instrument  resolution, 
making  it  impossible  to  verify  this  hypothesis.  Consequently,  we  chose  to  model  two  extreme  cases  for  the 
for  the  P:GeH  anisotropy,  A  =  pz/px.  The  first  assumes  no  change  in  p'z  compared  to  the  vacuum 
measurements  ( A  »  1  )  so  that  the  water  has  no  effect  on  the  OOP  resistivity.  The  second  case  assumes  a 
large  reduction  in  p'z  due  to  the  water  activation  (A  =  1). 

For  the  strong  anisotropy  (A  »  1)  case,  p'z  =  5  x  1014  Q  •  cm,  the  value  calculated  from  the  OOP 
measurements.  Here,  p'  =  2.3  x  ]  O'  Q  •  cm,  the  same  as  the  previously  modeled  smooth  sample.  The  result 
for  this  scenario  is  shown  in  figure  5(c).  Again  >99%  of  the  current  flows  through  the  topmost  layer.  The 
final  case  (figure  5(d))  again  assumes  ten  layers  of  water  activated  GeH,  but  this  time  the  water  activation 
is  taken  to  increase  the  interlayer  conduction  yielding  isotropic  resistivity  values,  A  =  1.  In  this  case  values 
of  Px  =  Pz  =  2.3  x  1  ()4  Q  •  cm  are  extracted  from  the  model.  Here,  >99%  of  the  current  travels  inside  of  the 
water  activated  region  and  is  equally  distributed  throughout  these  surface  layers. 

A  summary  of  the  modeling  results  is  shown  in  table  1 .  In  all  three  cases,  we  find  that  the  value  of  the 
in-plane  resistivity  is  bounded  2.3  x  103  <  p'x  <  2.3  x  104  Q  •  cm,  which  is  10  -  11  orders  of  magnitude 
smaller  than  the  calculated  out-of -plane  resistivity  value  of  pz  =  5  x  1 0 1 4  O  •  cm.  Additionally,  >99%  of 
the  current  is  confined  to  only  the  water  activated  layers  of  the  sample  in  all  of  the  modeled  cases. 
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Table  1.  Summary  of  modeling  results  for  the  cases  involving  water  activated  dopants. 


Number  of  Activated 
Layers 

Pz  (H  ■  cm) 

Px  (U  •  cm) 

A  —  Pz/ Px 

1 

5  x  1014 

2.3  x  103 

2.2  x  1011 

10 

5  x  1014 

2.3  x  103 

2.2  x  1011 

10 

2.3  x  104 

2.3  x  104 

1 

4.  Conclusion 

We  have  grown  both  doped  and  undoped  GeH  crystals  and  evaluated  their  transport  properties  for  the  first 
time.  By  incorporating  phosphorous  into  the  synthesis  of  the  CaGea  precursor,  it  is  possible  to  electronically 
dope  GeH,  leading  to  a  reduction  in  the  sample  resistance  by  up  to  three  orders  of  magnitude.  Electrical 
transport  measurements  conducted  on  bulk  crystalline  GeH  samples  using  two  different  contact  geometries 
show  a  large  resistance  anisotropy  between  the  in-plane  and  out-of-plane  directions.  The  resistivity  of  the 
crystals  exhibit  a  strong  atmospheric  dependence  with  the  resistivity  decreasing  by  up  to  three  orders  of 
magnitude  when  samples  are  exposed  to  air.  This  atmospheric  dependence  of  the  resistance,  along  with  the 
photocurrent  measurements  suggest  that  the  presence  of  water  vapor  in  the  air  is  activating  dopants  on  the 
surface  of  the  P:GeH  crystals.  To  better  understand  the  current  distribution  inside  these  layered  materials, 
we  developed  a  finite  element  model  which  indicates  that  nearly  all  of  the  current  travels  in  the  topmost 
water  activated  layers.  These  experimental  and  modeling  results  suggest  the  possibility  of  the  ability  to 
measure  exotic  single  layer  properties  in  the  top  layer  of  a  multilayered  group  IV  graphane  crystal  by 
making  electronic  contact  to  the  top  surface. 
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